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We study transport properties of a disordered tigiit-binding model (XX spin chain) in the pres- 
ence of dephasing. Focusing on diffusive behavior in the thermodynamic limit at high energies, 
we analytically derive the dependence of conductivity on dephasing and disorder strengths. As a 
function of dephasing, conductivity exhibits a single maximum at the optimal dephasing strength. 
The scaling of the position of this maximum with disorder strength is different for small and large 
disorder. In addition, we study periodic disorder for which we find a resonance phenomenon, with 
conductivity having two maxima as a function of dephasing strength. If disorder is nonzero only at 
a random fraction of all sites, conductivity is approximately the same as in the case of a disorder 
on all sites but with a rescaled disorder strength. 



I. INTRODUCTION 

Understanding quantum transport in simple sys- 
tems is of obvious relevance for understanding nano 
and mesoscopic devices as well as transport in real 
materials. One such model system is a tight-binding 
chain in the presence of disorder and dephasing, that 
has been studied already more than 30 years agoi. 
The goal is to understand the interplay of disor- 
der, incoherent processes and possibly interactions, 
on the transport properties. If the system is co- 
herent procedure is, at least in principle, straight- 
forward - one has to calculate the transmission by 
one of various approaches, see e.g. book @ for an 
overview. Non-coherent processes due to interaction 
with external degrees of freedom, for instance due to 
electron-phonon scattering, are more difficult to ac- 
count for. Correspondingly, they are also less under- 
stood. In the present work we study the influence of 
the environment, in particular its dephasing effects, 
on the bulk quantum transport in a disordered sys- 
tem. 

There are several approaches how to study in- 
coherent transport. A rigorous one in terms of 
nonequilibrium Green's function is difficult to an- 
alytically evaluate for all but the simplest sys- 
tems, while its numerical calculation is very time- 
consuming and limited to small systems. With that 
in mind, alternative, more phenomenological ap- 
proaches to incoherent transport are actively inves- 
tigated. One of the earliest suggestions accounts 
for the dephasing by introducing fictitious reser- 
voirs that break phase coherence^. This so-called 
Biittiker probe has been applied to the disordered 
tight-binding model* and is successfully used for the 
calculation of transmission through molecules, see 
e.g. Refs. d and H. A phenomenological model 
of dephasing can also be constructed by using an 
appropriate self-energy in a nonequilibrium Green's 



function^. Another approach consists of treating the 
system as being composed of statistically distributed 
coherent and incoherent parts (having an imaginary 
self-energy)^"^". 

In this paper we use different approach by writ- 
ing quantum master equation of the Lindblad type^— 
describing time evolution of the reduced density ma- 
trix of a tight-binding model without environment. 
Environmental dephasing is described in an effective 
way by dissipative operators that cause the decay 
of off-diagonal coherences (i.e., also of the current). 
Coupling the system in addition to reservoirs a true 
nonequilibrium steady state (NESS), reached after 
long time, is studied. Description with the Lind- 
blad equation has its advantages and disadvantages. 
Compared to other approaches, once one finds the 
NESS one has access to all many-body observables, 
i.e. to a complete state. Of course, this advan- 
tage can be harnessed only if one is able to solve for 
the NESS. Fortunately, a tight-binding model with 
dephasing is solvable'^^ in its Lindblad formulation. 
While a study of NESSs of Lindblad equations is 
a well-defined and important mathematical problem 
it has its down sides also. A number of approxi- 
mations, like weak coupling and fast-decaying bath 
correlations, is needed in standard derivations of the 
Lindblad equationi^. We note, however, that Lind- 
blad master equation is expected to provide good de- 
scription on timescales that are much larger than the 
timescale of bath correlations. It namely provides 
the most general description for quantum processes 
that preserve positivity and trace and form a dy- 
namical semigroup. Memory effects, not accounted 
for by the Lindblad equation, are expected to be 
smaller at higher temperatures, which is the regime 
that we study. In addition, we are interested in 
the bulk transport properties in the thermodynamic 
limit for which finite-size boundary effects are ex- 
pected to play no role. Compared to small systems, 



where more complicated master equations might be 
necessary for the correct description of transport, 
see e.g. Il4l . we can use simpler Lindblad equation to 
calculate bulk conductivity at high temperature. 

Tight-binding model with dephasing and disorder 
has recently become of interest also in the context of 
excitation transfer in biomolecules, in particular in 
photosynthetic complexes. There coupling with the 
environment, i.e., dephasing, increases the efficiency 
of excitation transfer because it counteracts local- 
ization due to disorder—"—. This, so-called environ- 
mentally assisted transport, can be in long homoge- 
neous chains explained by analytical results that we 
present. 

Our study proceeds in two steps. First, exploit- 
ing solvability, we transform the problem from the 
state space of density matrices, that is exponentially 
large in the chain length, to the one of two-point 
correlation functions, whose size is only quadratic in 
the chain length. We manage to analytically solve 
the resulting equations for some parameter regimes, 
while for other we resort to numerical solution. In 
addition to explaining the dependence of conductiv- 
ity on disorder and dephasing strength we also iden- 
tify an interesting resonance phenomenon for peri- 
odically placed disorder. 



II. XX CHAIN WITH DEPHASING 

We shall first present the model in the framework 
of Pauli matrices and show later a well-known fact 
that the model is equivalent to a nonequilibrium 
tight-binding model with an on-site disorder and ex- 
periencing dephasing due to environment. 

The Hamiltonian of the XX spin chain in an in- 
homogeneous magnetic field is given by 

L-l L 

j=i i=i 

with standard Pauli matrices and L being the chain 
length. We use units in which h = 1 and nearest- 
neighbor hopping is J = 1. We shall describe evo- 
lution of the system in the presence of environment 
by the Lindblad equation^ ^, 

^p^i[p,H] + /:^-(p)^£(p), (2) 

in which dissipative influence of the environment is 
described by the so-called dissipator C^'^. Action of 
a dissipator on an arbitrary operator p is expressed 



in terms of Lindblad operators Lk as, 

k 

Dissipator will in our case consist of two parts, 
£dis ^ £bath _^ £dcph^ rp^g g^g^ ^bath describes the 

action of two baths, one at each chain end, and will 
induce a nonequilibrium situation with current flow- 
ing through the system. The bath is written as a 
sum of a part acting only at the left end (site in- 
dex j = 1 and label "L") and a part acting only 
at the right end (site index j = L and label "R"), 
£bath ^ £bath £bath^ j,^^ Lindblad Operators at 

the left end are 

ll = Vr(i + M) <7t, 4 = v/r(i-/i) af, (4) 

while at the right end we have 

Lf=^T{l-ti) <j+, L« = yf(TT^(T-, (5) 

= {aj ± io'J)/2. Bath therefore flips the bound- 
ary spin up or down with certain probability. Pa- 
rameter /i plays the role of the forcing, trying to in- 
duce magnetization -t-/i at the left end and — /i at the 
right endi^. The other dissipative part, C^^p^^ repre- 
sents the influence of environment at each site, being 
for instance due to scattering on phonons. The de- 
phasing part £dcph _ ^'j'^^^ is a sum of 
each of which acts only at the j-th site and is de- 
scribed by a single Lindblad operator, 

Lf = (6) 

Due to dissipation the system's density matrix p{t) 
will after long time converge to a time-independent 
state Poo = limt^oo p{t) that is called a nonequi- 
librium steady state (NESS). For our model the 
NESS Poo is unique and can be formally expressed 
via the solution of Lindblad equation ([2|) as poo = 
limt_>oo exp (>Ct)p(0). In this work we shall be in- 
terested in the NESSs bulk transport properties in 
the presence of disorder and dephasing. Conductiv- 
ity is expected to be insensitive to details of bath 
Lindblad operators in the thermodynamic limit. All 
expectation values reported in the paper are with 
respect to the NESS state, e.g. (A) = tr (Apoo). 

Because we are interested in transport properties 
the central object we shall consider is the spin (i.e., 
magnetization) current jr, defined at site r through 
a continuity relation, resulting in expression jV = 
i [Ki ^r^r+i + giving 

> = 2(aM+i-«+i)- (7) 
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Dephasing alone causes an exponential decay of 
ofF-diagonal matrix elements (in the standard ba- 
sis in which cf" is diagonal). Therefore, due to 
£dephj'j^^ _ j£ |;]2ere would be no H and 

no £^'^**^ in the Lindblad equation, dephasing would 
cause an exponential decay of current with time, 
jrii) — exp (— 47i)jV(0). In the presence of Hamilto- 
nian and driving the evolution is more complicated, 
resulting in a nontrivial NESS. Evaluating the action 
(HI) of dephasing ([B]) on the identity and ct^ operators, 
we see that C^'''P\al) = £'i°P'^(lr) = hold. 

III. STATIONARY SOLUTION 

To find the NESS poo we have to find stationary 
solution of Eq. ([2]), i.e., solve £(/9oo) — 0. Expanding 
Poo in an operator basis the stationary equation can 
be written as a set of coupled linear equations for 
unknown expansion coefficients. Because the num- 
ber of unknown coefficients is 4^ — 1, and therefore 
grows exponentially with the system size L, it is for 
a generic system impossible to analytically obtain 
Poo. XX model with dephasing though has a nice 
property that these exponentially many equations 
split into smaller sets of equations that are uncou- 
ple d^^'^^1^^'^^ . Namely, there is a hierarchy of ob- 
servables according to the number of fermionic op- 
erators they contain. For instance, all two-point 
observables decouple from the rest, meaning that 
one can write a closed set of equations for observ- 
ables involving two fermionic operators. Once these 
are known, they can be used as a "source" term in 
equations for three-point observables, and so on for 
higher order correlations. For general consideration 
when such a hierarch y ap pears see Ref. [28 and 30. 

As shown in Ref. [12, one consequence of such 
structure is that the NESS can be written as 

L L+l-r 

^-E'E^^r<^ (8) 

r=2 ]=1 

where for r > 2 we define h''^^^^ = cr^ Z^^j^^^ + 

crJzj^7^''o-J+,,, with zf^ = a^j ■ ■ ■ a^^^^-^ being a 
string of r consecutive cr^'s. i/j^'' is defined as 
i/j^' = —(Jj- Let us stress that Eq.® is not just 



a perturbative expansion in p as it might seem at a 
first glance - it is an exact ansatai^i^ holding for any 
driving p. Or, in other words, terms proportional 
to p~ (and of higher order) not written explicitly in 
Eq. ([8]) are all orthogonal (using Hilbert-Schmidt in- 
ner product {A,B) = tr{A^B)) to operators in H 
and B. Expectation value of these observables is 
therefore obtained exactly to all orders in p by cal- 
culating only the expansion coefficients contained in 

(2) 

H and B. In particular, noting that ffr is a hop- 

(2) 

ping contribution to the energy density, Br current 
and Hr^^ magnetization, all quantities necessary to 
study transport are contained in the ansatz ([5]). 

Writing now £{poo) — using the ansatz ([H]), we 
get the aforementioned closed set of linear equations 

(r) (r) 

for unknown coefficients and bj . These are 

r + r/i^^^ - fef ^ = 0, 
-r + r4^'-6i'li = o, 

6f-6gi=0, j-2,...,L-2, (9) 
and for r > 2 

ih^-''~'^') + ihr''-'^')+ (10) 

+(6,-e,+._o/^r+Tf?^r=o 

+ (6,-e,+._i)fer-T^r^r-0' 

where T^-"^^ = + TSj^i + r6j+r-i,L- There are 
equations for exactly as many unknown coefficients. 
Solving Eas. pII|) therefore gives exact expectation 

(r) (r) 

values of observables Hj ' and Bj ' . 

The set of equations (|9ll0p can be compactly writ- 
ten in a matrix form by defining a hermitian corre- 
lation matrix C of size L x L with matrix elements 
C,,. ^ /^f-^+i) -Hifef-^+i) for k > J, ^ hf\ 
while Cj,k — C^^j for j > k. Let us define also 
^ = ii? — i J + TR, where i? is a diagonal disorder 
matrix, Ejj = e^, while nonzero matrix elements of 
J are Ji,i+i = 1, Ji+i.i = 1, while R has only two 
nonzero elements, Ri^i — Rn,n = 1- AH equa- 
tions (|9ll0p can now be written as a single matrix 
equation, 

AC + CA'' +2-fC ^ P, (11) 

where Pi i = —2T, Pn,n — 2r, while all other ele- 
ments of P are zero, and C — C ~ diag(C) is the 
correlation matrix without the diagonal. Physically, 
P represents the driving, E the disorder, R the cou- 
pling to baths, C term the dephasing while J is due 
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to XX Hamiltonian. Note that without dephasing, 
7 = 0, matrix equation (|lip would be of the Lya- 
punov type. Lyapunov equations appear in NESS 
solutions of quadratic systems, see e.g. [Sll 

Equation (|lip is the starting point for our study. 
In some cases we are able to solve it analytically, in 
other we resort to numerically exact solution using 
standard linear algebra packages that enable solu- 
tion for chain lengths L of upto several thousand. 

Expectation values of all operators and B - 
are trivially proportional to /i, see Eq. and we 
therefore without loss of generality from now on set 

= 1. Coupling strength to reservoirs T in general 
influences boundary resistance. Because we are in- 
terested in the regime of diffusive bulk transport, 
for which the boundary resistance does not mat- 
ter in the thermodynamic limit, we in addition set 
r = 1. For transport properties special attention 
shall be paid to the diagonal elements of the cor- 
relation matrix C, giving the magnetization pro- 
file, and to the imaginary part of the first near- 
diagonal, , giving the current. Note that cur- 
rent is due to Eq. © independent of the site index; 
henceforth we shall frequently use a shorter notation 

(2) 

b = bj \ resulting in spin current expectation value 
(j) = 46. Spin conductivity, being a proportionality 
constant between the gradient of a driving field, in 
our case magnetization, and a current, is therefore 
K = limi^oo / "^^^''L — lim^^oo 2i6. Here we used 
that in the thermodynamic limit (crf^ — erf ) — > 2 (for 
/i = 1). We see that the conductivity is finite and 
nonzero, i.e., we have a diffusive transport, provided 
the current scales as j ~ 1/L. 

Let us finally briefiy discuss the energy range to 
which the studied NESS states correspond. We pre- 
fer to use energy instead of temperature because 
thermalization in integrable systems, of which the 
XX chain is an instance, can depend on the choice of 
Lindblad reservoirs^-. There are two contributions 
to the energy ([T]), hopping and magnetization. As 
we shall see, the average magnetization profile is lin- 
ear. Each individual term ejcrf will have fluctuations 
of order ~ afi centered around zero, meaning that 
their sum gives negligible contribution to the energy 
density in the thermodynamic limit. Similar conclu- 
sion holds for ajcrj^i + crjcrj_|_;^. Namely, such term 

(2) 

is essentially given by ~ bj '{ej — ej+i)/j (see e.g. 
Eq. IA3[) . again giving negligible contribution to the 
energy density in the thermodynamic limit. NESS 
solution ([5]) is in the thermodynamic limit therefore 
always close to the identity matrix (i.e., an infinite 
temperature state). This means that we are study- 



ing transport properties at high energies. 

IV. FERMIONIC PICTURE 

It is helpful to translate quantities from spin 
language of Pauli matrices to the fermionic lan- 
guage of creation and annihilation Cj oper- 
ators, satisfying anticommutators {cj,Ck} = 0, 
{c|,cj,} = 0, and {cj,cj.} = <5j,fcl. Using Jordan- 
Wigner transformational^ Cj — — (o"! • • • cr|_i)cr^, 
Cj = -(ctj • • • cr|„;^)CT~, or its inverse, aj = 
-(af • • • +c]), aj = -iial ■ ■ ■ a|_i)(c, -c]), 

and a J = CjC^- — CjCj — 1 — 2nj, where we denote by 
Uj = CjCj a number operator at site j, the Hamil- 
tonian ([I} can be rewritten as a nearest-neighbor 
hopping, 

^ = E - c,4+i) + - 2^1.). (12) 

3 i 

Random magnetic fields translate to a diago- 
nal on-site disorder. Spin current translates to 
jr = 4i(cj!cr+i -I- Crc|!+i), i.e., is proportional to the 
particle current. Bath Lindblad operators (|H) on the 
other hand inject or absorb a particle, thereby in- 
ducing nonzero particle current through the system. 
Operators ijj'^'' and B^'p in the ansatz translate 
for r > 2 to 

i?f+^)=2i(c]c,+.+c,4+J. (13) 

They are two-point operators involving two 
fermions. Coefficients in the correlation matrix C 
PT|) therefore give expectation values of all two- 
point functions in the fermionic picture. In partic- 
ular, density profile is {tij ) = (1 + hf^)/2, while 

the particle (charge) current defined as = 

2i[nk,c\ck+i - CfcC^j] = i?^^-*, has expectation 
value, 

J ^ (jf = 26. (14) 
Conductivity in fermionic picture is therefore 

K = lim , '^■'k I ^ lij^^ i2b, (15) 

L^oo (n^ — ni) L->co 

and is the same as in spin language. In most of fig- 
ures we shall either show the dependence of particle 
current J^"\ simply denoted by j, or of conductivity 
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V. HOMOGENEOUS DISORDER 

In this section we study the simplest case in which 
the disorder is present at each site and has the same 
variance. Each is an independent random variable 
with a uniform distribution p{w) = const, in the in- 
terval w e [— \/3ct, \/3cr]. We therefore have disorder 
averaged values (ej) — and (e|) = a^. Our results 
presented do not depend on details of the distribu- 
tion p{w), only on its width a. Also, a homogeneous 
nonzero average (e-,) would not change any of the 
results presented, see discussion in Ref. 

In the next subsection we shall first study the case 
without dephasing for which one has Anderson local- 
ization (in the master equation setting). This regime 
is studied just to set the relevant length scale of lo- 
calization. New results about the interplay of de- 
phasing and disorder are then presented in the sub- 
sequent subsection IV Bl 



A. Localized phase, 7 = 

Without dephasing, 7 = 0, we have an ordi- 
nary tight-binding model with diagonal disorder in 
a nonequilibrium setting. One expects that An- 
derson localization^ of the conservative model, for 
which in ID all eigenstates of H are localized, will 
be reflected in an exponential decay of current with 
length L (at constant driving). This has been nu- 
merically observed in [271: here we determine the de- 
cay rate, i.e., the localization length, because it is 
an important length scale relevant also for 7 7^ 0. 
In Fig. [1] we can see that the localization length de- 
termined from the dependence of current scales as 
I ~ consistent with previous similar observa- 

tions in a nonequilibrium setting^. It is known from 
early studies of the Anderson model^^ that the dis- 
tribution of current is log- normal, i.e., distribution 
of its logarithm is normally distributed, away from 
band edges. For log-normal distribution the average 
current can be expressed as (j) = exp (m -I- s^/2), 
where m = (lnj),s = cinj- The question of con- 
ductance distribution in the Anderson model at the 
band center or at the band edge however is a rather 
delicate one. Namely, an otherwise universal single- 
parameter scaling theorySS^ saying essentially that 
the distribution function of conductance is a func- 
tion of a single parameter (say of its first moment), 
is modified in the band center— 1^'^'^^ . We can see 
in Fig. [5] that in our nonequilibrium setting the 
distribution of current is almost log-normal. The 
single-parameter scaling relation between the first 
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FIG. 1. An ensemble averaged logarithm of the particle 
current for different chain lengths L — 64, 128, 256, 512 
(bottom to top) and disorder strengths a, all for 7 = 0. 
We have a universal exponential decay due to localiza- 
tion, Inj ~ {—L/l), with localization length I ~ 
(solid line above numerical curves). Non-asymptotic be- 
havior for a^L < 10 is due to localization length being 
larger than L. 
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FIG. 2. (a) Distribution of particle current (points) is 
close to log-normal (full curve). All is for L = 128, a — 1. 
(b) Dependence of the variance of Inj on its average. 
Both plots are for zero dephasing, 7 = 0. 



two moments of Inj on the other hand does not hold. 
Whether small discrepancies visible in the decay of 
current in Fig. [T] are due to the band-center anomaly 
is not clear. 



B. Diffusive regime for nonzero dephasing 

We proceed to the main part of our work, that 
is to the case of nonzero dephasing 7. Relevant pa- 
rameters are disorder strength a, dephasing strength 
7 and the chain length L. Fixing a and 7 > 0, 
we find that the current always scales diffusive as 
j ^ l/L for sufficiently long chains. That is, in 
the thermodynamic limit a nonzero dephasing al- 
ways causes the system to become diffusive, regard- 
less whether it was ballistic (for cr = 0) or localized 
(for cr ^ 0) without dephasing. We note that in- 
teraction between particles still preserves diffusive 
nature of transport^!. We shall always discuss only 
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properties in the thermodynamic hmit, that is in the 
regime when the system is diffusive. For finite size 
effects that appear in the transition region from a 
locahzed/ballistic phase for small L's to a diffusive 
phase see Appendix 

We first focus on the case of large disorder a. For 
large cr the size of |Cjj+r| exponentially decreases 
with the distance r from the diagonal (we will see 
that this holds also for large 7). In the extreme case 
we can therefore approximate C with a tridiagonal 
matrix: we assume that the only nonzero elements 
of C are on the diagonal and in the 1st off-diagonal. 
Doing this approximation Eq. ljlip can be solved ex- 
actly for any L, a and 7, see Appendix [X] In the 
thermodynamic limit the expression for h is 
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L 272 + a2 ' 
immediately giving conductivity (fT5 

1 



7 + 



(16) 



(17) 
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Transport properties in condensed-matter systems 
are many times discussed in terms of scattering 
lengths. Heuristically, one often assumes that dif- 
ferent scattering processes are independent and that 
one can simply sum- up scattering rates, i.e., recipro- 
cal scattering lengths. Interpreting analytical result 
for conductivity (|17|) in such a way one can say that 
the total scattering length ^tot has an independent 
contribution from dephasing ^dcph = 1/7 and disor- 
der Zdis = 27/ct2, giving 1/Ztot = l/^h + l/'dis with 
K oc Itot- This of course holds only in the regime 
of validity of Eq. p^ . that is for large cr or large 
7. In Fig. [3] we can see that formula p7|) agrees 
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FIG. 3. (Color online) Dependence of conductivity k, on 
dephasing strength 7 for large disorder a. Numerical re- 
sults (full red curve) agree well with theoretical formula, 
Ea. (fT7)l (dotted curves). 

well with numerical results for large tr. Two gross 



features are visible: (i) ^ I/7 dependence of k for 
large dephasing ~ as the dephasing becomes large, 
scattering due to impurities can be neglected, and 
(ii) for small dephasing the transport is dominated 
by a break-down of localization caused by nonzero 
dephasing; correspondingly, k increases with the de- 
phasing. A combined effect of the two contributions 
causes a maximum in conductivity at an interme- 
diate dephasing strength. Such behavior has been 
observed in other studiesis^"— of dephasing effects 
on transport in a disordered tight-binding model 
as well as in transport properties of molecules or 
molecular aggregate o^i^^'^^i^^ . The location of the 
maximum in k scales according to (jl7p as 7niax = 
cr/v^, with the value of k at the maximum being 
'«max = l/(cr\/2) = l/(27inax)- Arouud (7 « 4 and 
smaller, appreciable differences between numerical 
results and the theoretical Eq. P3|) are visible. The 
reason for this failure at smaller a is that the ap- 
proximation of the correlation matrix by a tridiag- 
onal matrix becomes increasingly worse. For even 
smaller a it completely fails. Numerical results in 
this regime of cr are shown in Fig. 21 As opposed to 
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FIG. 4. Conductivity at smaller disorder. Large-a the- 
ory (dotted curves), Eq. I|17p . still works for sufficiently 
large dephasing 7 ^ cr. 

large cr, the location of the maximum in the regime 
of small cr scales as 7max ^ cr^. Fig. m while the con- 
ductivity at the maximum is Kmax ^ l/o"^ ^ l/7max- 
We see that the scaling of 7max changes with cr. 
This can be nicely seen in Fig. [5] The crossover be- 
tween the two asymptotic behaviors happens around 
cr « 3. We observe that this coincides with the dis- 
order strength at which localization length for 7 = 
becomes / w 1. 

Observe in Fig. |4] also that formula ([T7| still holds 
for sufficiently large 7, i.e., for dephasing larger than 
about ss cr. For 7 < cr and small cr, where (|17l) does 
not work anymore, we were not able to obtain the- 
oretical prediction for the current and k. Including 
more that three diagonals in the calculation does not 
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FIG. 5. Location of the maximum 7max in i^('y) as a 
function of disorder strength a. Lines are asymptotic 
theoretical scahng, having two regimes, points are nu- 
merics. 



bring a significant improvement on Eq. (|17p . It seems 
that the problem is in this regime strongly domi- 
nated by fluctuations and all elements of C have to 
be taken into account. Perturbation theory in ei- 
ther 7 or a is also not successful, see discussion at 
the end of Appendix [Cj In Fig. [5] we show in more 
detail this small-cr behavior. We can see that the de- 
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FIG. 6. (Color online) Scaled conductivity k for small 
disorder a. 

pendence on 7 is rather complicated to the left of the 
maximum at 7max- Overall, k scales as ~ 1/^^ ^^'^ 
is approximately a function of the scaled variable 
j/cr'^- Scaling though is not perfect (see the region 
left of the maximum), perhaps pointing to a nontriv- 
ial (non-additive) interplay of dephasing and disor- 
der. Whether this complicated behavior for small 
a is in any way related to anomalous properties of 
the Anderson model at the band cente r^^'^^ (i.e., at 
zero energy) for weak uncorrelated disorder remains 
to be explored. Our Lindblad reservoirs namely in- 
duce a NESS with energy close to 0, i.e., in the band 
center. 



VI. DILUTED DISORDER 

In this section we discuss the case when disorder 
is absent at some sites. Such situation of a diluted 
disorder is relevant for real materials that can be 
prepared with different concentration of impurities. 
On theoretical side, it is known that in a Hamilto- 
nian system non-homogeneous disorder, for instance 
a correlated one, can have a profound effect on the 
Anderson localization. In ID correlations in disor- 
der can namely cause a delocalization, see e.g. H^. 
It is also worth mentioning that the band center 
can be special in terms of localization properties for 
weak (uncorrelated) diagonal2^i^ as well as for off- 
diagonal disorder—. 

We shall study two cases: (i) randomly placed di- 
luted disorder in which a fraction p of randomly cho- 
sen sites have disorder, and (ii) periodically placed 
disorder for which disorder is placed at every A-th 
site. We will see that the transport properties of 
case (i) are very similar to the case of a homogeneous 
disorder with a renormalized disorder strength ^^pa. 
Case (ii) though is qualitatively different. Disorder 
distribution function shall in all cases be the same 
unform distribution as before. 

Before going to the properties of nonequilibrium 
states we briefly comment on the localization prop- 
erties of a corresponding conservative system. Cal- 
culating eigenvalues and eigenstates of a tridiagonal 
matrix describing a tight-binding model with disor- 
der at randomly placed fraction p of sites, case (i), 
all eigenstates are still localized. Their localization 
length though differs. A fraction p of eigenstates 
has localization length that is the same as if disor- 
der would be at every site, while (1 — p)L eigen- 
states have larger localization length and are clus- 
tered around the band center at E — 0. For periodic 
disorder, case (ii), delocalized states appear at the 
band center—. We found that all eigenstates with 
eigenenergies \E\ < (for A = 2) are delocalized. 
We expect that this fraction ~ 1/L of all eigenstates 
will strongly influence transport properties out of 
equilibrium if dephasing is zero, 7 = 0. 



A. Randomly placed diluted disorder 

From the 3-diagonal derivation in Appendix [XJ 
holding for large 7, we see that the fact that dis- 
order is not present at each site will be reflected in 
an additional prefactor in front of a^. Instead of 
La^ as for the homogeneous case we will have pLa^, 
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while there is no disorder at ah other sites. 



FIG. 7. (Color online) Conductivity for a diluted dis- 
order. Disorder is placed on random pL sites. Black 
dotted curves are for a homogeneous disorder (p = 1) 
with strength g^i, blue and red are for a diluted case 
with o" = a^f/^/p. 



resuhing in the conductivity 



1 



(18) 



The equation is the same as would be for the ho- 
mogeneous case with an effective disorder strength 
fcf = v^f- This expression is expected to be valid 
for large 7. For smaller dephasing, even though we 
do not have a theoretical formula, we can still try to 
approximate the diluted situation with a homoge- 
neous one at the effective disorder strength inferred 
from large 7 behavior (jl8p . In Fig.[7]we compare re- 
sults for a diluted disorder with the ones for a homo- 
geneous case with a smaller effective disorder. The 
agreement between k for a randomly diluted disor- 
der and a homogeneous one with a rescaled effective 
disorder is not perfect to the left of the maximum. 

That for diluted disorder the disorder strength 
is just renormalized is not very surprising. Such 
result follows from the assumption that scatter- 
ing events are independent and therefore scattering 
length scales as ^ 1/p. Scaling of such form has 
been observed experimentally in ID spin chain ma- 
terials described by the isotropic Heisenberg model 
in which scattering event are due to impurities^S. 
We expect that findings of the present manuscript 
for the XX chain would not qualitatively change in 
the presence of interaction, i.e., for the Heisenberg 
model. 



Periodic disorder 




L=32, 64, 128 



FIG. 8. Periodic disorder with period A = 2 and 
7 = 0. Scaling of the particle current j with o for 
L = 32, 64, 128. Difference between (a) and (b) is only in 
the scaling on the a:-axis. (a) For La^ < 10 (when local- 
ization length / is larger than L) periodicity is irrelevant 
and the dependence is the same as for a homogeneous 
disorder with a rescaled disorder strength a/\/2 at every 
site (black solid curve "hom."; the same exponentially 
localized data as in Fig. [l|. Scaling variable is a^L. (b) 
For larger disorder the scaling variable is gL and the 
decay is algebraic, j ~ l/(Lcr)'^ with a ^2. 



First, we check conductivity for 7 = 0. Results 
of solving Ea. (jlip and averaging over disorder are 
shown in Fig. [8] For sufficiently small disorder cr, 
such that the localization length / is larger than L, 
the dependence of current j on a is the same as for a 
homogeneous disorder with a renormalized disorder 
strength of size a^Jl/X. One can see in Fig. [8^ that 
the differences between an exponentially localized 
homogeneous case and the periodic case begin to 
appear only for a^L ^ 10. For larger disorder, when 
periodicity becomes important, the nature of decay 
with L changes (Fig. Current begins to decays 
in an algebraic way as j ~ 1/ L^. Such decay can be 
traced back to properties of the Hamiltonian system, 
which has a fraction 1/L of delocalized eigenstates 
having eigenenergies of order 1/L located around 
the band center. 
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In this subsection we place disorder of strength a 
at every A-th site, i.e., at sites j = kX — 1, fc S N, 



FIG. 9. Current for disorder on every 2nd site and the- 
oretical formula (1211). dotted curves. L — 256. 
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We next move to nonzero dephasing for which the 
transport is always diffusive in the thermodynamic 
limit. System size at which this diffusive behavior is 
reached is larger than in the homogeneous case and 
depends on o and 7. Main difference with respect 
to the homogeneous or randomly diluted case is that 
the dependence of ^(7) has two peaks. An additional 
resonant peak appears at smaller 7 for sufficiently 
large disorder cr, see Fig. [9] While the right peak 
is correctly described by Eq. dTS]) with p = 1/A, we 
have not been able to theoretically describe the left 
peak. Empirically we find that its position is 



71 « l/(3Aa), 

while its height scales as ^ l/c^ 
position of the right peak is from 



(19) 



Note that the 



72 = cr/v 2A 



(20) 



(while its height is ^ l/f) and therefore scales differ- 
ently than the left peak. The left peak is essentially 
a resonant phenomenon while the right one is due to 
competition between disorder and dephasing. 

Even though we are not able to theoretically de- 
scribe the left peak we can derive the location of 
the minimum between the two peaks. Assuming 
that nonzero elements of C are only on the diag- 
onal and two off-diagonals {b^^\h^'^^ and b'^^\h^^'>), 
we can solve for b, see Appendix [Bj The result is in 
the thermodynamic limit 



1 



^7(1+2^) 



(21) 



where cjcf — a/y/X. Note that the difference from 
the 3-diagonal approximation ()f 6|) is only in f in 
the denominator. One can see in Fig.[S]that Ea. (pT|) 
correctly describes the minimum. Lo cation o f the 
minimum is at 47^;^^ = cr^f — 2 — dof^/ cr^j — 8; this 
minimum is real for act > V& ~ 2.83 while its lo- 
cation goes toward 7min ^ ■\/3/2 ~ 0.73 for large a. 
We therefore expect the left resonant maximum at 
71 to appear only for ct > 3 which is indeed con- 
firmed by numerical calculation shown in Fig. 1101 
From Eq. (PT|) we also see that the height at the 
minimum scales as x l/cr^ which is the same as the 
empirical scaling of the left maximum. To also de- 
scribe the maximum, and not just the minimum, we 
would have to solve for at least the diagonal and 4 
off-diagonals in C, which proves to be too difficult. 

In Fig. [TU] we can see that locations of the two 
maxima is nicely predicted by Eas. ([TO)) and (PO)) . 
Dependence of k for periodic disorder, having two 




0.0001 0.001 



FIG. 10. (Color online) Density plot of a scaled conduc- 
tivity a K as a function of 7 and a for disorder on every 
A = 2 site (L = 256). We multiply conductivity by a in 
order to have better contrast for two peaks. Two black 
curves denote position of the maxima, left curve Ea. (ll9|l . 
and right curve Eq. (|2Up . 




0.0001 0.001 



FIG. 11. (Color online) Scaled conductivity an as a func- 
tion of 7 and a for disorder at randomly chosen half of 
sites (p — 0.5, L — 256). Black curve is theoretical 
prediction for the location of the ma^dmum at large a, 
Eg. pH) with A = 2. 



maxima, can be compared to the one for randomly 
placed disorder (with the same average density) , for 
which only the right maximum is present, see Fig. 1111 



VII. CONCLUSION 

We have studied conductivity in a tight-binding 
model with on-site disorder and experiencing de- 
phasing. Transport is induced by a nonequilibrium 
driving described by the Lindblad master equation. 
Writing linear system of equations for all two-point 
fermionic expectation values in the steady state we 
study dependence of conductivity on dephasing and 
disorder strength. For large disorder or large de- 
phasing we derive an analytic expression showing 
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that the conductivity has a single maximum as a 
function of dephasing. We also discuss the case 
of diluted disorder. For randomly placed diluted 
disorder conductivity is approximately the same as 
for a homogeneous disorder with a rescaled disorder 
strength. For periodically placed disorder though a 
second resonant maximum in conductivity appears, 
exhibiting different scaling on parameters as the first 
maximum. We acknowledge support by the grant 
Pl-0044. 



Appendix A: 3-diagonal approximation 

For large disorder o (and any 7) or for large 7 
the size of correlations |Cj_j-|_r| rapidly decreases 
away from the diagonal, that is, one has |Cjj+r| ^ 
exp {—r/p) with small p. In such a case one can ap- 
proximate C by its diagonal Cj,j and its two nearest 
diagonals, Cj j+i, Cj+i.^ . Setting all other matrix 
elements of C to zero one can solve Eq. ([TT|) exactly 
as follows. 

We use standard notation for the unknown terms 



in C, h 



(1) 



im(a 



and h 



(2) 



Re(Cjj+i). Two nontrivial equations (|lip on the 
diagonal are 



2Th^ 
2Vh\ 



(1) 



2b = -2T 
2b = 2T, 



(Al) 



while off-diagonal we have imaginary parts, 



bk+{h[^^ 
267+ (4'' 



£2 







4") + 4'^(£2-e3) = 



267 + (/if ^ 



6A:+(/iW,-/i^^0 + 4-i(^L-i-^L) = 



(A2) 



,(1) 



(2) 



where A: = F -I- 27, and real part, 

- kh^^^ + 6(ei - €2) 

.(2 



2^hY> +b{e,-e,+i)^0 



khf_^ + b{eL-i - el) = 0. 



(A3) 



To get the coefficient of the current b we sum all 
equations (|A2[) . and use the relation /ii^"* - 
2(5 - r)/r obtained from This gives 



,(1) 



5(2fc+(L-3)27) + ^ hf{e,-e,+,) + -{b-V) = 0. 
i=i 

(A4) 

It is worth pointing out that Eq. (|A4[) holds true in 
general, not just in the case of a 3-diagonal approx- 
imation. Inserting into this relation expressions for 



obtained from Eqs. (|A3p . we get 

27 



c + 272(L- 1) 



e.+i)72 



(A5) 



with the boundary term c that is independent 
of the length L being c = 27(F + i) + 
[(ei - 62)2 + (ei_i- 6^)2]. 
Let Ej be distributed according to a distribution 
having a finite 2nd moment, (e^) = a^. Then 
ej — Ej+i has 2nd moment equal to 2a^. Accord- 
ing to the central limit theorem a sum of terms 
(ej — Ej+i)' will converge to a Gaussian distributed 
random variable with a nonzero average 2(L — 3)a^ 

and a variance oc L. Expression Ej=2^(^j ~ ^i-i-i)'/2 
therefore becomes increasingly sharply peaked about 
its average, with relative fluctuations being of order 
^ l/\/Z. Therefore, in the limit L — 00 when we 
can neglect boundary terms, we can write 



27 



L(272 + ^2) 



(A6) 



Appendix B: 5-diagonal approximation for 
periodic disorder 

Let us have disorder only on odd sites, e2j = 0, 
and assume that the correlation matrix is 5-diagonal, 
Cj,j+r = for r > 2. The goal would be to go be- 
yond the 3-diagonal approximation in order to de- 
scribe the 2nd resonance peak visible for instance in 
Fig. [9l Although we will fail at that - one would 
need to solve at least a 11-diagonal approximation 
to get the 2nd peak (which we are though not able 
to do analytically) - we will nevertheless correctly 
describe the minimum between the two peaks giving 
us at least some insight. 

We will neglect boundary effects and solve equa- 
tions in the leading order in 7/172. ^pj-^jg jg justified 
because the 2nd resonance maximum appears only 
for large a. First, we note that because e2j = 



10 



one 



has 



0. Then we also observe that 

there- 



2j+i 7/0" and can be neglected. We 
fore start with the fact that ft,'-^-' are zero on even 
sites while b^^^ are (approximately) zero on odd sites. 



On even sites we have b, 



(3) 



-2, - -ih^S - /^^r-/i)/(27). 
Writing now equations for the real part of the near- 
diagonal at two consecutive sites in the bulk, we have 



,(2j) 



'2^ 



,(2) 



27 



27 



e2j+ib = 



£2^+16 = 0. 



From these two equations we see that only two 
nearest-neighbor h^'^^ are coupled. Instead of L cou- 
pled equations we have a set of uncoupled 2x2 equa- 
tions. Solving them we get 



For instance, in Fig. [T^ we can see that at 7 = 
10~^ one needs L = 1024 or larger, while at 7 = 10 
chain length L = 32 is enough to observe diffu- 
sive behavior. The system size L at which diffu- 
sive behavior is reached is smaller for larger disorder 
strength a, as can be seen in Fig. [T31 

We also note that the perturbation theory in 7 
around 7 = has a convergence radius that shrinks 
to zero with the system size L. The failure of a naive 




,(2) 



^2j 



27 



(27)2 



-,<^2j+ib. 



(Bl) 



Inserting these into Eq. (jA4p we get in the thermo- 
dynamic limit L — >■ 00, 



FIG. 12. Scaled particle current for a — 0.5 and dif- 
ferent sizes L. Convergence to diffusive scaling, seen as 
an overlap of curves, is for smaller 7 reached only for 
large sizes (data for L = 1024 and 2048 overlap and are 
indistinguishable on the scale of the plot). 



27 



L(272 + ct2 



(B2) 



1+1/(272). 



This expression correctly describes the minimum be- 
tween the two peaks in Fig. [HI although it fails to re- 
produce the left maximum. Nevertheless, the scaling 
of the conductivity at the minimum is the same as 
the empirically found scaling at the left maximum. 



Appendix C: Finite-size effects 

Limits L — )■ 00 and 7 — > do not commute. De- 
pending on the order of the two we get either a diffu- 
sive behavior j ~ 1/L or an insulating j ^ due 
to localization. We are interested in the diffusive be- 
havior obtained in the correct thermodynamic limit 
of letting L — >■ 00 first. In order to observe diffusive 
behavior one therefore has to take sufficiently large 
L for each particular 7. 




FIG. 13. Same as Fig. [12 but for a = A. Data for 
L = 128, 256, 512 overlap. 



perturbation theory to predict diffusive behavior at 
small nonzero dephasing is not very surprising be- 
cause the nature of transport changes abruptly at 
7 = 0. Similar are problems trying to apply per- 
turbation theory in a in order to obtain behavior at 
small a and 7 (for instance in Fig. |6l). 
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